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Abstract 

Interest to the tasks by bound with systems of the quasilinear 
equations recently has revived. Therefore we have decided to expose 
in electronic sort our operation published in Translated from Teo- 
reticheskaya i Matematicheskaya Fizika, Vol. 101, No. 2, pp. 189-199, 
November, 1994. Original article submitted October 26, 1993. 

Exact single-wave and two-wave solutions of systems of equations 
of the Newell- Whitehead type are presented. The Painleve test and 
calculations in the spirit of Hirota are used to construct these solu- 
tions. 



INTRODUCTION 

In [1-4], methods were developed for constructing the explicit expressions 
that describe nonlinear waves (kinks) that are solutions of the semilinear 
parabolic equations 

U t -U XX = F(U), (0.1) 

where F{U) are quadratic and cubic polynomials. The solutions of Eq. (0.1) 
considered in the cited studies satisfy boundary conditions that are deter- 
mined by the roots of the nonlinearity, i.e., solutions propagating in strips 
were constructed. Besides bounded solutions, solutions having a disconti- 
nuity of the second kind (so-called monsters) were considered. In addotion, 
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expression were obtained for solutions that describe the interaction of kinks 
and also kinks in discontinuous solutions. At the present time, only two 
types of solution of Eq. (0.1) represented by explicit expressions are known 
- so-called self-similar (single-phase, single-wave) solutions and solutions de- 
scribing the interaction of self-similar waves (two-phase solutions) [1,4]. The 
self-similar solutions have the form (for the notation, see below) 



'V'exp(Ar) 

For the nonself-similar solutions describing the interaction of single-phase 
waves, there is a representation of the form 

V ? (exp(Air 1 ),exp(A 2 r 2 ),ri,r 2 ) 

X{n,r 2 ) = U{— - - -), 

^(exp (AiTi), exp (A 2 r 2 ), t u t 2 ) 

(0.3) 

Here,U(z) is a function that is entire in some region Q e C 1 ; ip and ip 
are polynomials; r = x + pt + cis the self-similar variable; A and p are 
phase (related) constants; and c is an arbitrary constant. Similarly, Tj = 
x + Pi + Ct, i = /, 2, where Aj and Pi - are phase constants, q is an arbitrary 
constant, and the phases ,Tj are independent, i.e., a linear combination of 
them with integer coefficients does not vanish identically. We note that self- 
similar solutions of (0.1) are determined by the solution of the ordinary 
differential equation 

\VX ~ AY - Fix) = 0- (0.4) 

It is shown in [3,4] that for equations of the form (0.4) the store of solutions 
of the form (0.2) and, hence, the set of constants A and p is very restricted. 
Similar results have also been obtained by means of the Painleve test [5-7], 
the idea of which is to represent the solution of an equation of the form (1) 
by a Laurent series in the neighborhood of a moving pole, i.e., to expand 
the required function in powers of r, the self-similar variable. The Painleve 
test, applied to various equations of the form (0.1), showed that the set of 
possible velocities for self-similar solutions represented by a Laurent series 
is, as a rule, small and can be uniquely determined. In addition, for each 
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concrete value of the velocity p obtained by means of the Painleve test the 
Laurent series can uniquely, up to a shift constant, determine a function of 
the form (0.2) . Therefore, in the case of a positive Painleve test, depending 
on its results, it is natural to assume the existence for equations of the form 
(0.1) of solutions in the class of functions defined by the expressions (0.2)- 
(0.3). In this case, one can substitute in the original equation the general 
form of the assumed solution (the ansatz), where <p and ip are polynomials 
with undetermined coefficients. Equating the coefficients of the same powers 
of the exponentials or (and) equal powers of the self-similar variables, we 
obtain are overdetermined, in general, system of nonlinear algebraic equa- 
tions for the unknown coefficients of the polynomials if and ip and the phase 
constants. The solution of this system (if it exists) is necessary and suffi- 
cient for constructing an explicit expression for the solution of the original 
equation. We note that it was precisely an assumption about the form of 
the solutions of equations of the form (0.1) that made possible the successful 
use of constructive methods of solution finding, for example, modification of 
Hirota's method ([2], p. 184,[2'],p.68). 

In this paper, we describe a modification of the Painleve test for sys- 
tems of semilinear parabolic equations. It turns on that the set of possible 
velocities of self-similar solutions of systems is, as in the case of individual 
equations, restricted and can be uniquely determined. The analogy with the 
scalar case made it possible to assume that the systems have solutions of 
the definite form, for the finding of which the method of undetermined co- 
efficients described above was used. It is characteristic that in the case of 
the systems considered in this paper, as in the one-dimensional case, it is, 
as a rule, possible to find all similar solutions whose velocity is determined 
by the Painleve test by making a restriction to functions of the form (0.2). 
By means of the same method we have succeeded in constructing two-phase 
solutions of the form (0.3) of some of the consider systems. 

We note however that the mechanism of interaction of waves that are 
solutions of systems of semilinear parabolic equation has been found to be 
very specific. Thus, in [4] (p. 44,[2']p.70) there is an asymptotic description 
of the process of wave creation from a finite perturbation for the Kolmogorov 
- Petrovskii - Piskunov (KPP) - Fisher equation, from which there follows 
the impossibility of constructing an expression that describes the interaction 
of the waves in the form dictated by Hirota's method. It turned out that the 
situation is quite different for systems of semilinear parabolic equations. In 
the paper, we consider example of a system for which explicit expressions for 
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a solution describing the annihilation of kinks are constructed (In this text 
the figures are not reduced. Look figures in the original of a paper in TMF 
. see Fig. 2 below). 

The subject of the investigation in the present paper is systems of the 
form 

u t -u xx = u\i-u m -e n ), 

e t - e xx = -Bu k e\ (0.5) 

where U, 9 are the unknown functions,!? is a nonvanishing constant, and 
m, n, k, q, I are natural numbers. The paper consist of two sections. In Sec. 
1, we give self-similar solutions of a system that is one of the forms of the 
Belousov - Zhabotinskii reaction model [8,9] and has, in addition, a biophys- 
ical application [10]. In Sec. 2, we give self-similar and two-phase solutions 
of various systems describing processes of the type of nonlinear kinetics. 

1 SELF-SIMILAR SOLUTIONS 

IN THE BELOUSOV-ZHABOTINSKII 
REACTION MODEL 

Suppose that in (0.5)m = n = k — q — I — l,B^0,B^ 1/2, i.e., the system 
has the form 

u t -u xx = u(i-u-e), 

9 t -9 xx = -BU9, (1.1) 
Then the following theorem holds. 

Theorem 1.1 All self-similar solutions of the system (1.1) that can be rep- 
resented by a Laurent series in a punctured neighborhood of the moving pole 
x = —pt + c, where p is the velocity and c is an arbitrary constant, have the 
form 

L r_ 1 e _ (l-g)(2expr + ex P 2r 

(1 + expr) 2 ' (1 + expr) 2 

T = ax + b t ^ a 2 = 5/6, b =-55/6; (1.2) 
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_ 1 - 2 exp r) (1-B) exp 2r 

(1 + expr) 2 ' (1+expr) 2 ' 

r = ax + bt, a 2 = -B/6, b = -55/6, p = b/a; (1.3) 

Proof. 

The fact that the given functions solve the system can be verified by 
substitution. We consider in more detal question of the uniqueness of the 
presented self-similar solutions associated with modification of the Painleve 
test for system. We shall seek a solution in the form of Laurent series in 
powers of r. 

oo 

U = ^2 a k+qi (x + pt + ci) fe , 

k=—qi 

oo 

9= B k+q2 (x + pt + Cl ) h , (1.4) 

k=-q 2 

where q\ and q 2 are certain natural numbers that determine the order of the 
poles of the functions U and 9, respectively, and in what follows we shall, 
without loss of generality, assume that the constants c\ and c 2 are equal to 
zero. In the general case, the possible values of q% and q 2 are determined 
by means of Newton's polygon method; however, in the present case it can 
be done by using elementary arguments and equating the exponents of the 
definitely highest negative powers of the second derivatives and the quadratic 
terms of the polynomials. It is easy to see that in the given case qi = q 2 = 2, 
i.e., both functions U and 9 have a pole of second order. Substituting the 
obtained series in the system and equating the coefficients of equal powers of 
r, we obtain a system of nonlinear algebraic equations for the determination 
of the coefficients of the series: 

a (a + b - 6) = 0, b (a B - 6) = 0, 
2aia + CLib — 2a x + a bi — 2a p = 0, 
a x b Q B + a bxB - 2b x - 2b p = 0, 

2a 2 a + a 2 &o + «i 2 + ai&i — &iP + ^0^2 — a = 0, 
a 2 b B + aibiB + a b 2 B — b±p — 0, 
2a 3 a + a 3 b + 2a 2 ai + a 2 b\ + a x b 2 — a\ + a bs = 0, 
B(a 3 b + a 2 bi + a x b 2 + a bs) = 0, 
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2a 4 a + CI4O0 — 2a 4 + 2a 3 ai + a 3 b\ + a 3 p + a\ 

+a 2 6 2 — a2 + ai& 3 + a &4 — 0, 
a 4 &o-B + 0301-8 + a 2 o 2 5 + aib^B + a »4-B - 26 4 + b 3 p = 0, 
2a 2 a + a 5 6 — 6a 5 + 2a 4 a x + a 4 6i + 2a 4 p + 2a 3 6 2 

+a 3 6 2 — a 3 + a 2 6 3 + ai6 4 + a b 5 = 0, 
a 5 b B + a A b\B + a 3 b 2 B + a 2 b 3 B + a\b 4 B + a o 5 5 — 6b 5 + 2b 4 p = 0, 
2a 6 a + a 6 6 — 12a 6 + 2a 5 ai + a 5 6i + 3a 5 p + 2a 4 a 2 

+a 4 6 2 - a4 + a 3 2 + a 3 6 3 + a 2 6 4 + ai6 5 + a o 6 = 0, 
a 6 b B + a^b\B + a 4 b 2 B + a 3 b 3 B + a 2 b 4 B + aib 5 B 

+a o 6 5 - 126 6 + 3b 5 p = 0. (1.5) 

Solving successively this system for the unknown coefficients a k and bk, 
we obtain 

6 6(5-1) 6p 6(5 -l)p 

a °-B> ~ 5 ' ai "55' 61 " 55 ' 
255 -p 2 (1 - 5)(p 2 + 255) 

° 2 - 505 ' 62 " 505 ' 



a 3 



p 3 . (5 - l)p 3 



2505 ' 2505 



1255 2 - 7p 4 (5- l)(1255 2 -7p 4 ) 

O4 — — , Oa — — , 

50005 ' 50005 

13755 2 - 79p 4 , p(l-5)(13755 2 -7V 
a 5 = ^tt^t; > h = 



750005 750005 
375006 6 - 6255 2 p 2 + 36p 4 

° 6 " 37500(5 - 1) • (L6) 

It is readily seen that in the final equation of this system the coefficient 
of the unknown b$ is equal to zero and the obtained equation is an expres- 
sion from which the possible values of the velocity can be determined. This 
expression is usually called the dispersion relation. It has the form 

(25-/)(36p 4 - 6255 2 ) = 0. (1.7) 

By virtue of the conditions of the theorem, the solution of 5g.(1.6H) for p 
gives two (up to the sign) possible values of the 

±5^ ±5iVB 

^ = ^g-. ^ = ~jr- (L8) 
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For each particular value of the velocity p , the coefficients and bk of the 
series for U and 9 are uniquely determined; therefore, there exist only two 
different self-similar solutions of the system that can be represented by a 
Laurent series with second-order pole. It is obvious that to each of these 
solutions there corresponds one of the two possible values of the velocity 
and that for a definite value of the parameter B with fixed sign these values 
are different. By virtue of the invariance of th substitution x— > —x in 
the system (1.5), reversal of the sign of the phase constant a corresponds 
to reflection of the traveling wave with respect to the ordinate with change 
of the direction of its motion, i.e., reversal of the velocity. Thus, the above 
solutions completely exhaust the store of velocities obtained by means of the 
Painleve test. □ 



Remark 1.1 It is clear that for any B ^ only one of the two solutions 
(1.2) and (1.3) is real. The boundary condition that real solutions of the 
system (1.1) satisfy have the form 

U x ^OO * 0, dx^OO * 1 B, Ux^ — OQ * 1; ^X^ — CO * 0. (l'^) 



Remark 1.2 It follows from Eq. (1.6) that for B = 1/2 there may also exist 
other solutions in addition to (1 .2) and (1 .3) that are given by explicit 
expressions. 

Note that as B — > 1 in each of the solutions (1.2) and (1.3) 6 — > 0, and 
at the same time the solution of the system it continuously transformed into 
the solution of the scalar equation 

Ut-U X x-U(l-U) = 0. (1.10) 

which is known as the KPP - Fisher equation. The Painleve test for Eq. 
(1.8) gives, up to the sign, two possible values of that velocity for solutions 
that can be represented by a Laurent series, namely 

Pl - 2 = 7i' P3 ' 4= vT (L11) 

The solutions with such values of speed were earlier calculated in [2] p. 184. 
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It is therefore obvious that there do not exist any other solutions of Eq. 
(1.8) that can be expanded in a Laurent series apart from those that are 
obtained from the expressions (1.2) and (1.3) for B= 1 . In addition, it is 
clear that in the given case the velocity (1.7) are continuous functions of the 
parameter B of the system on the complete domain of definition, i.e., for any 
B^O will be shown below that this is not always the case. 

2 SELF-SIMILAR AND TWO-PHASE SO- 
LUTIONS IN MODELS OF NONLINEAR 
KINETICS 

Suppose that in (0.5) m — n — k — 2 : q — I — 1. We obtain the system 

u t -u xx = u(i-u 2 -e 2 ), 

9 t -9 xx = -BU 2 9, (2.1) 

As in the case of the system (1.1), the order of the poles of the functions U 
and 9 in the system (2.1) is determined. In the given case, ql = q2 = 1 in 
(1.4) and (1.5). The unique (up to sign) possible expression for the velocity 
that is obtained by me; of the Painleve test has the form 

It proved to be possible to construct solutions of the system (2. 1) for B = 2. 

Theorem 2.1 All self-similar solutions of the system (2. 1), B = 2 , repre- 
sented by a Laurent series in the punctured neighborhood of the moving pole 
x = —pt + c , where p is the velocity, and c is an arbitrary constant, have 
the form 

U= ±X , 9 ±6Xpr 



1 + exp r ' 1 + exp r ' 

r = ax + bt, a — ±1, b = — 1, p = b/a = ±l; (2.3) 

U = ^-, 9 ±6Xpr 



a — t 1 + exp r ' 

T = ax + bt, b = 2a, p = b/a = 2; (2.4) 
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where c is an arbitrary constant. Proof. The fact that functions (2.3) and 
(2.4) satisfy the system can be verified by substitution. The Painleve test for 
the considered system gives two possible values of the velocity for solutions 
that can be represented by a Laurent series: pi = 2,p 2 = ±1. It is with 
such velocities that the above solutions move. From this the statement of 
the theorem follows, and the theorem is proved. □ 

Remark 2.1 The value of the expression (2.2) for B = 2(p = y/3) is different 
from the values of the velocity given in Theorem 2 for the solutions of the 
system (2.1) for the same value of the parameter B. Thus, the expression 
(2.2) is not a continuous function of the parameter B of the system (2.1). 

Proposition 1. The system (2.1), B = 2, has a two-phase (nonself- 
similar) solution of the form 

ai+expr 2 ±n 

U = , V = , 

0,1 + ti + exp r 2 ai + n + exp r 2 

Ti = dix + bit, t 2 = a 2 x + b 2 t, 

bi = -2ai, b 2 = 1, a 2 = 1; (2.5) 

where ai is an arbitrary constant. The proposition is verified by substitution. 
The fnteraction of the self-similar waves (2.3) and (2.4) described by (2.5) is 
analogous to the interaction of a kink with a discontinuous solution, which 
can be described by explicit expressions for individual equations with cubic 
nonlinearity [4] (p. 39-42), [2']p. 51. Figure 1 gives the graph of the function U 
for ai = 1. It can be seen that on the collision of the singularity with the front 
of the kink the order of the latter increases. For en = — 1 , the interaction of 
the discontinuous solution with the kink described by the function U leads 
to annihilation of the singularity and the formation of one kink. A situation 
analogous to the considered cases also holds for the function 6. We turn to 
the next example. Suppose that in (5) m — n — q — I, I — k — 2. The 
resulting system has the form 

U t -U xx = U 2 (l-U-6), 

6 t - 9 XX = -BU 2 9, (2.6) 

In this case, it proved possible to construct both self-similar and two-phase 
solutions containing the parameter B of the original system (2.6) as a free 
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parameter. The Painleve test for the system (2.6) gives two possible values 
of the velocity for self-similar solutions that can be represented by a Laurent 
series, namely 

Pi = ^2 p 2 = ±V2B. (2.7) 

Note that this is true for any B ^ 1/3, since for B = 1/3 the dispersion 
relation vanishes identically. A consequence of this result is the following 
theorem, the proof of which is completely analogous to that of Theorems 1 
and 2. 

Theorem 2.2 All self-similar solutions of the system (2.6), B ^ 1/3, that 
can be represented by a Laurent series in the punctured neighborhood of the 
moving polex = —pt + c, where p is the velocity and c is an arbitrary constant, 
have the form 

a (1- B)a 2 expr 

U = , t> = , 

ai + a 2 exp r C\ + a 2 ex P r 

T = ax + bt, a 2 = 5/2, b=-B/2, p = b/a = ±^L, (2.8) 

v2 

ci,C2 are arbitrary constants; 

u= e _ VB(l-B)r 
rVB + aV^' rVB + aV^' 
T = ax + bt, b= -aV2B, p = -V2B, (2.9) 

a is an arbitrary constant. 

Remark 2.2 As in the case (1.1), for 5 = 1/3 there can be other solutions 
in addition to (2.8) and (2.9) that can be represented by explicit expressions. 

Proposition 2. The system (2.6) has a two-phase solution of the form 
u _ ai + a 2 exp r 2 ^ _ a 2 (l - B)n 



ai + a 2 (r 1 +expr 2 ' a x + a 2 (ri + exp r 2 ' 
ri = a±x + bit, t 2 = a 2 x + b 2 t, 

6i = -1a x a 2 , a\ = B/2, b 2 = B/2, (2.10) 
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is ai arbitrary constant. The validity of the proposition is proved by sub- 
stituting the expressions (2.10) in the system (2.6). As in the case of (1.1), 
the self-similar solutions of the system (2.6), and also the two-phase solu- 
tion (2.10) describing their interaction are reduced continuously in the limit 
B — > 1, respectively, to the self-similar and two-phase solutions of the one- 
dimensional problem known as one of the forms of the Zel'dovich equation 
([2], p. 73, p. 198): U t -U xx -U 2 {l-U) =0. Here, as in the system (1.1), the 
velocities (2.7) are continuous functions of the parameter B on the complete 
domain of definition. Note that all the given solutions of the system (2.6) are 
real for B > 0. The functions (2.10) describe the interaction of self-similar 
solutions of the system (2.6). The nature of this interaction is similar to the 
case considered above (see Fig. 1 in the original of a paper in TMF ). We 
now give an example of the system (5) in which it proved possible to find not 
only self-similar solutions but also a bounded two-phase solution that is the 
result of interaction of waves propagating in strips. 

Suppose that in (5) m — k — 2, q — n — I — I. Then we obtain the 
system 

u t -u xx = u(i-u 2 -e), 

9 t -9 xx = -BU 2 9, (2.11) 



Theorem 2.3 The system (2.11) has nontrivial self-similar solutions that 
can be represented by a Laurent series for a unique value of the parameter 
B, namely, B = 1. 

Proof. Functions U and 6 that can be represented by a Laurent series 
in powers of t and satisfy the system (2.11) have a first-order pole, i.e., 
the coefficients of a and bo must be nonzero. It follows from the system of 
equations for these coefficients that this is possible only if B — 1: 

a (a 2 -2)=0, b (a 2 B -2) = 0. (2.12) 

The theorem is proved. □ 

Proposition 3. There exist self-similar solutions of the system (2.11), 
B — 1, having the form 

U= ±X , 9= 6Xpri , 
1 + exp r 1 + exp r 2 

r = ax + bt, b=-l/2, a 2 = 1/2, (2.13) 
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Proposition 4. There exists a two-phase solution of the system (2.11), 
B — 1, having the form 

1-expn expr 2 

^ = T~i 1 ' y = T7 1 ' 

1 + exp Ti + exp r 2 1 + cxp r x + exp r 2 

ri = + bit, t 2 = a 2 x + b 2 t, 

bi = 0, ai = V2, a 2 = 1/V2, b 2 = -1/2, (2.14) 

□ 

The validity of both propositions can be proved by substituting the given 
solutions in the system (2.11). The two-phase functions U and B are the re- 
sult of interaction of the self-similar waves described in Proposition 3. Graphs 
of the two-phase solution are shown in Fig. 2. It can be clearly seen that 
the interaction described by the function 6 (see Fig. 2a) represents the an- 
nihilation of kinks propagating in the strip between and 1 (for t > to, the 
solution is close to zero). In the case of the function ?7(see Fig. 2b), there 
is an interaction of the type of merging of waves propagating in different 
strips [for t > to, a single wave close to the stationary solution-of Semenov's 
equation ([2], p. 193) is formed]. 
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